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Abstract 

We investigate radiative corrections to K — > 3tt decays. In particular, we extend the 
non-relativistic framework developed recently to include real and virtual photons and 
show that, in a well-defined power counting scheme, the results reproduce corrections 
obtained in the relativistic calculation. Real photons are included exactly, beyond the 
soft-photon approximation, and we compare the result with the latter. The singularities 
generated by pionium near threshold are investigated, and a region is identified where 
standard perturbation theory in the fine structure constant a may be applied. We 
expect that the formulae provided allow one to extract S-wave irir scattering lengths 
from the cusp effect in these decays with high precision. 
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1 Introduction 

The investigation of the so-called cusp effect in K — > 37r decays [1-3] has become 
a fully competitive method for the extraction of S-wave nn scattering lengths 
from experimental data. The theoretical basis for the cusp analysis was provided 
in Refs. [3-7] (see also Ref. [8]), whereas experimental results were reported in 
Refs. [9-13]. It turns out that the accuracy of the experimental value for the dif- 
ference ao — a2 of 7T7T scattering lengths 1 is mainly limited by two shortcomings in 
present theoretical descriptions [3-7] of the decay amplitudes: missing contribu- 
tions from n > 3 loops, and missing radiative corrections. This article is devoted 
to an evaluation of the latter. Leaving aside for the moment the contributions 
from higher loops, we expect that K — > 3ir decays, with radiative corrections in- 
cluded, and combined with the information gained from K e 4 decays [14-17] and 
the pionium lifetime [18, 19], have the potential to test the very precise theoretical 
prediction of the scattering lengths [20,21]. 

We briefly describe the framework in which we perform the calculation. Gener- 
ically, we provide algebraic expressions for the decay spectra, including the effects 
of real and virtual photons, in a coherent approximation scheme which respects 
the strictures of unitarity and analyticity. The pertinent expressions for the decay 
spectra contain several free parameters, to be adjusted such that the experimen- 
tal distributions are reproduced. Two of these parameters are the S-wave nn 
scattering lengths ao,2 we are after. 

The cusp effect in both K + — > 7r 7r°7r + and Kl — > 3tt° is seen in the decay 
spectrum with respect to the invariant mass squared S3 of 7r°7r° pairs. In order 
to obtain an infrared-finite quantity, one has to take the radiation of an addi- 
tional final-state photon into account as soon as electromagnetic corrections are 

1 We use the notation ai for the (dimensionless) S-wave scattering lengths with isospin i". 
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included, 



dT_ 

ds 3 



dT(K -> 3tt) rfr(K -> 3tt 7 ) 



+ 0(a 2 ). (1) 



Here, a = 1/137.036 denotes the fine structure constant, whereas the maximal 
energy of photons not resolved explicitly in the experiment is denoted by E max 
[measured in the kaon rest frame]. The mentioned approximation scheme is 
achieved by performing the calculations of the decay spectra in the framework 
of non-relativistic effective field theory (EFT) [6,7], extended in this article to 
include photons. Non-relativistic EFT produces a correlated expansion in im 
scattering lengths a and in the non-relativistic momentum parameter e, chosen 
such that the pion 3- momenta in the final state are counted according to |p| = 
(9(e). This expansion is supplemented here with a further expansion parameter 
e 2 = 4na. 

In the following, we will present the photonic corrections to the decay spectra 
up-to-and-including terms of order e 2 a°e 4 , e 2 a 1 e 2 for the two "main modes" Kl — > 
3n° and K + — > 7r°7r°7r + that show a cusp behavior in the invariant 7r°7r° mass 
spectrum inside the physical region. On the other hand, in the two "auxiliary 
modes" K L — > n + n^n and K + — > 7r + 7r + 7r~, we restrict the accuracy up-to-and- 
including terms of order e 2 a°e 4 . Furthermore, we confine ourselves to the so-called 
"soft-photon approximation" in the fully charged channel. [The framework is set 
up in such a manner that the pion masses are not affected by contributions from 
virtual photons, as a result of which it is fully consistent to set all pion masses to 
their physical values from the very beginning. No expansion in e 2 is performed 
in the pion mass difference.] 

Our results can then be represented in the form 



dT 

ds 3 



dT int 

= n(s 3 ,E m ^)— (2) 

E 7 <E max ^ 



for the three channels involving at least one neutral pion in the final state, and 
dT\ = n++-( Sl , s 2 , s 3 , E max )dT int (3) 

for K + — > 7r + 7r + 7r~, where both the decay spectrum dF int /ds 3 and the differential 
decay width dT mt are given in terms of a squared matrix element |.M mt | 2 . Here, 
fl(s 3 , -Emax) denotes a channel- dependent correction factor which is due to real and 
virtual photons hooked to charged external legs. All other photonic corrections 
are collected in the modified matrix element Ai mt . We will provide both, the 
correction factor and Ai mt , for all four decay channels. 

Our work is not the first one to consider radiative corrections to K — > 3ir 
decays - recent examples include Refs. [22-27]. Whereas Isidori's article [25] fits 
into the framework considered here, the calculations of Refs. [22-24] are of a wider 
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scope and not useful to extract tvtt scattering lengths with high precision by use 
of the approach proposed in Ref. [3], a framework also adapted here. We refer 
the interested reader to Sect. 9 for a more detailed comparison of our calculation 
with these works. 

The article is composed as follows. We start out with a brief recapitulation 
of the non-relativistic EFT framework without photons in Sect. 2. The inclusion 
of photons on the basis of non-relativistic effective Lagrangians is discussed in 
Sect. 3. Section 4 treats the essential photon loop diagrams for the analysis at 
hand, showing how they can be calculated in non-relativistic EFT, from which we 
derive general power counting rules in Sect. 5. In Sect. 6 we discuss the special 
role of pionium for the cusp region. Results for the correction factors fl(s 3 , -E max ) 
and amplitudes Ai mt are given in Sect. 7. This section contains the main results 
of our investigation. We discuss the various sources of theoretical uncertainties 
for the determination of the mr scattering lengths a — a 2 in Sect. 8. In Sect. 9, we 
compare our work to other articles dealing with radiative corrections in K — > 3n 
decays. Finally, we close with a summary. Some background material is relegated 
to the appendices: Appendix A collects our notation and explains the necessary 
kinematics, Appendices B-D contain the calculation of radiative corrections in 
a relativistic theory. Appendix E is a sample calculation clarifying the interplay 
between ultraviolet and infrared divergences in the threshold expansion, and a 
discussion of the specific symmetry properties of the Kl — > 3n° phase space can 
be found in Appendix F. 

2 Non-relativistic effective theory in the ab- 
sence of photons 

We now describe the covariant non-relativistic framework that will be invoked 
later on. In this section, we recall the method in the case where electromagnetic 
corrections are discarded [6, 7]. The inclusion of photons will be discussed in the 
following sections. 

2.1 Counting rules 

The complete Lagrangian of the effective theory is Ck + £wn, where Ck contains 
K + — > 37r, K L — > 37r vertices, and C wn describes elastic nn scattering in the final 
state. We omit 6-pion couplings, see later. The interactions are ordered according 
to specific counting rules [6,7]. A formal small parameter e is introduced, and 
the various quantities are counted as follows. 

• The masses Mr, Mk l , M n , M^o are counted as 0(1). 

• All three-momenta (spatial derivatives) are counted as 0(e). 

• Kinetic energies are counted as C(e 2 ). 
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As a result of this, one has to further count 

• M K - 3M n = 0(e 2 ), M Kl - 3M n = 0(e 2 ). 

• — M n o =D(e 2 ). 

Here, M K) M Kl ,M w , M w o denote the masses of K + , K L , tt + , 7r°, respectively. 

We shall treat the pion-pion interaction perturbatively, owing to the smallness 
of the pion-pion scattering lengths. For bookkeeping reasons, we introduce an 
additional formal small parameter a, which stands for the pion-pion scattering 
length. All four-pion couplings are regarded as quantities of order a. At the 
end, the amplitudes are given in the form of a double expansion in e and a. 
This expansion is correlated, since adding one more pion loop with an additional 
four-pion vertex also increases the power of e by one. 

In the following, we provide the Lagrangians necessary to carry out the cal- 
culation of the amplitudes for K + — > 3ir, K L — > 3n. 

2.2 The Lagrangians 

We start with the mr interaction and consider the following five physical channels 
in ttV* -> n c n d : (aft; cd) = (1) (00; 00), (2) (+0; +0), (3) (+-; 00), (4) (+-; +-), 
(5) (++;++). [We omit the channel it~it — > it^it , because this amplitude is 
identical to 7r + 7r° — > 7r + 7r° by charge invariance.] The Lagrangian takes the form 

5 

£™ = 2£ $IW± (idt - W±) $± + 2®IW Q (idt - Wo) $ + Yl A > ( 4 ) 

± i=i 

where $ m , m — ±, are non-relativistic pion field operators, and W m = 
a/ M 2 m — A, with A the Laplacian. At leading order in the non-relativistic ex- 
pansion, one simply finds 

C i = x i C i (<Z>i& d <Z> a <Z> b + h.c.) + ... , (5) 

where the ellipsis stands for terms of order e 2 and higher. Higher-order terms 
are explicitly given in Refs. [6,7]. The low-energy constants Cj are matched to 
the physical threshold amplitudes below in Sect. 3.2. To simplify the resulting 
expressions, we have furthermore introduced the combinatorial factors x\ = = 
1/4, x 2 = X3 = Xi = 1. Finally, we note that we omit local 6-pion couplings. 
Their contribution to the K — > 3n amplitudes is purely imaginary in the non- 
relativistic framework, and of order e 4 . 

Next, we consider the non-relativistic Lagrangians that describe K — > 3n 
decays in different channels. The following notation is useful: 

= {V%tf , {V% = (W K , <V) , (Kl), = {V%K* L , (P% = (W K , <V) , 

(6) 
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for n = a, b, c, d, where K,K L denote the non-relativistic fields for the K + ,K L 
mesons, and W K = \J M\ — A, W K = M] <l — A. The Lagrangians that 
describe K + — > 3ir and K L — > 3n decays at tree level are given by 

C K = 2K^W K (id t - W K )K + X - Go (K*$+$l + h.c.) 



+ 2 



+ - H„ (K*<t>_$l + h.c.) 



+ i ffl {((^:_M^_)^ +A .,} + ..., ( 7 ) 



t K = 2K\W K (iQ t - W K )K L + L„(K[ *„* + 4>_ + ft.c.) 

Ma 



+ ^„(^ + k)+..., (8) 

where the ellipsis stands for the higher-order terms in e. The couplings Gj, H i: L i; 
fTj are assumed to be real. 

2.3 Loops 

The Lagrangians displayed above generate pion loop contributions, according to 
the standard Feynman diagrammatic technique. The pion propagator is given by 

S Up) = 1 , \ / x 1 o Tf ' ( 9 ) 

2w m (p) w m (p) -p° -lO 



where w m (p) = a/ M 2 ™ + p 2 . The charged and neutral kaon propagators Sk(p), 
Sk(p) are defined similarly. 

However, since the heavy scales (pion and kaon masses) are explicitly present 
in the particle propagators, the simple and straightforward counting rules in e, 
which were introduced at tree level, will be ruined by loop corrections. As it is 
well known (see, e.g., Refs. [28-30]), the counting rules can be restored, imposing 
additional prescriptions (the threshold expansion) built on top of the Feynman 
rules. In our case, these prescriptions state: 

• For a given Feynman integral, expand all inhomogeneous factors (e.g. square 
roots) in the propagators in powers of the inverse pion masses. 

• Integrate the expanded Feynman integral term by term. 
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Figure 1: Two topologically distinct non-relativistic two- loop graphs describing 
the final-state nir re-scattering in the decay K — > 3tt, with = (pi +^2)^- 

• Re-sum the final result to all orders of the inverse pion masses. 

It should be noted that the framework is more complicated technically (due to 
the presence of the square roots that ensure the correct relativistic dispersion 
law for the particles in the loops) than the standard non-relativistic EFT (see, 
e.g., Ref. [30]). However, the bonus is that the location of the low-energy sin- 
gularities in different amplitudes is Lorentz-covariant. This is an important ad- 
vantage in the study of three-particle final states, where different two-particle 
sub-systems are in general not in the center-of-mass frame. Note also that, with- 
out re-summing the final result, we again arrive at the standard formulation of 
the non-relativistic EFT. 

The procedure outlined above has been carried out up to two loops [6, 7, 31]. 
Details of the calculations will be provided in Ref. [31]. We only give a short 
summary here. 

The one-loop contributions are proportional to the basic integral 
f d D l ?'A 1/2 fs M 2 M?) 

■Us) = / 4^ s a (p - i)s b (i) = lgifrA> , (io) 

where s = P 2 , and X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx denotes the triangle 
function. We see that, in contrast to the conventional non-relativistic framework, 
the basic integral is explicitly covariant. 2 

There are two types of two-loop diagrams with 4-pion vertices, shown in Fig. 1. 
At the order we are working, it suffices to consider the diagrams in Fig. 1 with 
non-derivative couplings only. In this case, the contributions of both diagrams 
depend only on the variable s, where 

s = Q\ Q^ipi+PiY- (11) 

The non-trivial contribution from Fig. la is proportional to 

M(s) = S wk 5 W? 5 Sa(P ' 1 ~ Q Sb ® s °W Sd ( Q - k) (12) 

2 This result can also be obtained in the framework developed in Refs. [32, 33]. 
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A short discussion of this integral is given in Ref. [6]. There, it is shown that one 
may write 

M(s) = F(M a ,M b ,M c ,M d ;s)+..., (13) 

where F is ultraviolet finite and contains the full non-analytic behavior of the 
two-loop diagram in the low-energy domain, whereas the ellipsis denotes terms 
that amount to a redefinition of the tree-level couplings in C K , C K and which are 
therefore dropped. A one-dimensional integral representation for F is provided 
in Ref. [6]. The relevant integrals can be also performed analytically [7]. 

Finally, the diagram in Fig. lb is a product of two one-loop diagrams given 
in Eq. (10). 

3 Including photons 

3.1 The Lagrangian with photons 

In this and the following sections we discuss the inclusion of photons in the 
non-relativistic framework. Since power counting is more complicated in the 
presence of photons, we first consider examples of simple vertices and diagrams 
with photons, without trying to order them according to the power counting. 
The kinetic part of the Lagrangian after minimal substitution takes the form 

£kin = -\ F (lv F >lv + ^(4AW±$i - i(AW±$ ± ) f $ ± - 2$| t Wl$ ± ) 

+ 2$JW (i3t - W )% + (iK*D t W K K - i{D t W K K)^K - 2K^W 2 K Kj 

+ 2K{W K (id t -W K )K L . (14) 
Here, F^ = d^A u — d u A^ denotes the electromagnetic field strength tensor and 
A$± = (dt T ieA )<Z>± , D t K = (d t - ieA )K , 

W± = y/Ml - D 2 , W K = \jM\ - D 2 , 

D$ ± = (V±ieA)$±, DK = (V + ieA)K . (15) 

In an explicit calculation of the tree-level photon-pion and photon-kaon vertices, 
the roots W±, Wk are understood to be expanded in the inverse of the corre- 
sponding mass. After resummation of all terms proportional to e, the relativistic 
result is explicitly reproduced at tree level, 

(^ip'Wm^ip)) = ±e(p'+p)», (K + (p')\j™(0)\K + (p)) = e(p'+p),. (16) 

jem ^ eno ^ es electromagnetic current, calculated from the Lagrangian in a 
standard manner. Note that the non-relativistic Lagrangian Eq. (14) does not 
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Figure 2: Effective-range coupling of the Coulomb photon to a charged pion. The 
coupling is proportional to the pion charge radius squared. 



contain a term for pion pair creation by a photon, 7 — > 7r + 7r _ , which would 
involve a high-energy photon. This omission eliminates all one-photon reducible 
graphs from the K — > 3tt amplitudes considered in the following. 

The spatial derivatives in the four-pion vertex are replaced by the covariant 
derivatives in the presence of photons. In addition, non-minimal terms, contain- 
ing the electric E and magnetic B fields, will be present. At lowest order, the 
non-local coupling of the Coulomb photon field A with e.g. the pion emerges, 
see Fig. 2. The pertinent coupling constant is proportional to the mean charge 
radius squared of the pion [30]. However, at the order we are working, this 
term in the Lagrangian can be eliminated by using the equations of motion and 
thus it does not affect any physical observable (we refer the interested reader to 
Ref. [30] for a detailed discussion of this issue). It can also be shown that higher- 
order non-minimal couplings do not contribute at the order we are working. In 
the following, we discard all non-minimal couplings in the Lagrangian from the 
beginning. Further, we also discard the contribution arising from vacuum po- 
larization generated by an electron-positron loop (see again Ref. [30]). We have 
checked that this correction is tiny and amounts to a few percent change of the 
(already small) Coulomb correction in the "auxiliary" decay modes. 

The loop calculations in the non-relativistic theory are most easily done in 
the Coulomb gauge. Here, first the time component Aq is removed from the 
Lagrangian by using the equations of motion. This procedure generates the non- 
local operator 

in the Lagrangian. In the following, for ease of understanding, we keep calling 
the pertinent diagrams as "generated by the exchange of Coulomb photons." The 
propagator of the transverse photons in the Coulomb gauge is given by 

/1 / P W \ 

d i xe ikx (0\TA i (x)A j (0)\0) = -— U ij - . (18) 

Below, the virtual photon corrections to the strong amplitudes are calculated at 
C(e 2 ) (with the exception of the discussion in Sect. 6), i.e. we include one virtual 
photon at most. In order to obey counting rules, as in the purely strong case, one 



10 



has to perform threshold expansions. In the following section we explain how this 
can be done on the basis of several examples. The counting rules in the presence 
of photons are discussed in full generality afterwards in Sect. 5. 

3.2 Matching 

The effective coupling constants in the Lagrangian can be related to physical 
observables in the underlying relativistic field theory. This procedure goes under 
the name of matching. Here, chiral perturbation theory is understood to be the 
underlying theory, which we bear in mind while performing the matching. 

In order to determine the (four-pion) coupling constants in the non-relativistic 
theory through matching, one has to calculate the pion-pion scattering amplitude 
both in the non-relativistic theory and in the underlying theory. These amplitudes 
are then expanded in the pion momenta at threshold, and the pertinent expansion 
coefficients are set equal (up to a given order). One has to distinguish the isospin 
symmetry limit and the case with broken isospin symmetry. 

3.2.1 Isospin symmetry limit 

In the isospin symmetry limit, the expansion of the relativistic nn scattering 
amplitude at threshold reads 

ReTi(s,t) = Ai + 0(e 2 ) . (19) 

The bar indicates the quantities in the isospin-symmetric world, in which the 
pion has the mass M n = 139.57 MeV. In terms of the standard scattering lengths 
a and a 2 , one has 

3^! = N(a + 2a 2 ) , 2A 2 = Na 2 , 3A 3 = N(a 2 - a ) , 

6A 4 = N{2a + a 2 ) , A 5 = Na 2 , N = 32n , (20) 

with a = 0.220 ± 0.005, a 2 = -0.0444 ± 0.0010, a - a 2 = 0.265 ± 0.004 [21]. 
Still in the isospin symmetry limit, the couplings Cj are related to these threshold 
parameters according to 

2Q = A, , (21) 

where we have dropped higher-order terms in the threshold parameters. 

We wish to note that the matching condition is universal in the isospin- 
symmetric case and determines the coupling constants in terms of the effective- 
range expansion parameters only. There is no explicit reference to the Lagrangian 
of the underlying relativistic theory. 
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3.2.2 Isospin-violating case 



The isospin-breaking corrections to the isospin-symmetric result emerge from 
electromagnetic effects as well as from the up/down quark mass difference and 
have the following general form (see, e.g., Refs. [29,30]): 

Ci = Ci + hi(m d - m u ) 2 + h 2 e 2 + C(e 4 , e 2 (m d - m u ) 2 , (m d - m u ) 4 ) , (22) 

where the coefficients h\ >2 depend on the quark mass m — \ (m u + m d ). Note that 
the corrections are no more universal, and in order to carry out the calculations, 
one has to use the explicit form of the underlying Lagrangian of the relativistic 
theory. Further, in order to perform the matching in the presence of photons, 
one has first to remove the infrared-singular piece of the amplitude at threshold. 
The procedure is described in Refs. [29, 30] and will not be repeated here. 
At leading order in chiral perturbation theory one finds [29,34] 

2671,2,5 = Ai,2 l5 (l - V), 2C 3 = A 3 {1 + 77/3), 2C 4 = A A {1 + 77), (23) 

where r] = (M 2 — M 2 )/M 2 = 6.5 x 10~ 2 . 

The calculations at one- loop level have been carried out, in particular, for the 
process 7r + 7r~ — > 7r°7r° that multiplies the cusp in the K — > 37r decays at lowest 
order. The final result can be extracted, e.g., from Eqs. (4.14), (4.28) and (4.29) 
of Ref. [29], 

327T 

2C 3 = — [a -a 2 + (0.61 ±0.16) x 10~ 2 ] 

= 4,(1 + (2.3 ±0.6) x 10~ 2 ) , (24) 

where in the second line we have used a — a 2 = 0.265. The bulk of the total 
correction is already given by the tree-level result in Eq. (23), and most of the 
uncertainty stems from the so-called "electromagnetic" low-energy constants in 
the 0(e 2 p 2 ) chiral Lagrangian [34]. The higher-order corrections to the matching 
of Ci, 2,4,5 are less important, because these couplings enter the cusp amplitude at 
the sub-leading order. In general, one may conclude that isospin-breaking correc- 
tions beyond tree level - given in Eq. (23) - amount to a systematic uncertainty 
in the Ci at the percent level. 

4 Examples of loop calculations with photons 

4.1 Ultraviolet and infrared divergences: general remarks 

The loop diagrams with photons in the non-relativistic EFT have both ultraviolet 
and infrared divergences. Their treatment is different: whereas the ultraviolet 
divergences are removed by renormalization, the infrared ones cancel in the total 
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decay rate, including the decay into soft photons. Further, the non-relativistic 
EFT exactly reproduces the structure of the infrared divergences in the underly- 
ing relativistic theory. On the other hand, the high-energy behavior in these two 
theories, in general, differs. 

Throughout this paper, we tame both ultraviolet and infrared divergences by 
dimensional regularization. In addition, the threshold expansion is used in the 
Feynman integrals in order to ensure the validity of the counting rules in the 
presence of photons. One has to face, however, the following problem. Both 
ultraviolet and infrared divergences show up as poles in the amplitudes in the 
variable D as D — > 4. The origin of the singularity (either ultraviolet or infrared) 
can be easily identified, and one may attach a subscript to these poles, writing 
{jjzi) uv an d (7534) IR explicitly. At the order we are working, it can be checked 
that, if threshold expansion is not applied, the infrared poles that emerge in the 
non-relativistic EFT, as expected, are identical to those in the pertinent graphs 
of the underlying relativistic theory. 

However, the threshold expansion generally changes the behavior of the dia- 
grams in the infrared, and the structure of the infrared poles does not match the 
underlying theory anymore. In order to see this, we note that applying thresh- 
old expansion and the subsequent use of the no-scale argument in dimensional 
regularization amounts to the replacement 
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D-4J UV \D-A /IR 



(25) 



in the pertinent terms. For more details, we refer the reader to Appendix E, 
where we illustrate the issue with the vertex diagram. 

In this paper, we adopt the following simple solution of the problem. As 
becomes clear from the above expression, the coefficient of the infrared pole 
changes by a low-energy "polynomial" as a result of the threshold expansion. For 
this reason, we do not differentiate between infrared and ultraviolet divergences 
(detaching the subscripts "UV" and "IR" everywhere) and remove all remaining 
divergences by the counterterms of the Lagrangian in a threshold-expanded EFT. 
This leads to the correct result for all observables, since all infrared divergences 
cancel at the end. 



4.2 Self-energy of the pion 

The self-energy diagram of the charged pion shown in Fig. 3a is given by the 
following integral 

Only the transverse photons contribute. 
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Figure 3: Examples of loop graphs with photons in the non-relativistic EFT. 
Double and single full lines denote kaons and pions, respectively, wiggly lines 
stand for Coulomb, curly lines for transverse photons, a: pion self-energy (trans- 
verse photon); b: vertex diagram with Coulomb photon; c: vertex diagram with 
transverse photon; d: vertex diagram in crossed channel (transverse photon); e: 
Coulomb photon exchange in the strong loop. 



In order to calculate S(p), we perform the Cauchy integral over 1° and further 
use the identity 

1 1 _ 1 1 

2w±(p - 1) w±{p -\)-p° + |1| ~2w±(p - 1) w±(p - 1) +p° - |1| 



T4(p-i)-(p»-|i|) a " 

Since both w±(p — 1), p° = 0(1) in the e counting, the threshold expansion in 
the first term generates polynomials in the integration momenta. Using no-scale 
arguments, it is seen that the integral of the first term vanishes in dimensional 
regularization. 

The denominator in the second term can be written as f2 + 2p°|l| — 2pl, where 
= tf±(p) — (p ) 2 = M 2 — p 2 = 0(e 2 ). In order to expand the denominator, 
one has to establish the power counting rules for the photon 3-momentum 1. One 
observes that in this particular diagram 1 should be assigned the power 0(e 2 ), 
otherwise one arrives at no-scale integrals after expansion in e, which vanish in 
dimensional regularization. Assuming 1 = 0(e 2 ) and expanding this denominator 
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in the last term, which is of order e 3 (the first two terms are of order e 2 ), we find 3 



e 2 p 2 f Q. y~ 2 / In 2tt^ 2 
2w±(p)p° V2p°y I 1 rfJ(27r) d r(! 



S(p) = ^rr^o(^n) ( 1 -l)7^7^ r ( rf - 1 ) r (2-rf) + ..., (28) 



where d — D — 1, and the ellipses stand for higher-order terms in the e-expansion. 
In the limit d — > 3 we obtain 



3w±(p)(p ) 2 V 16vr 2 V /up 3 



(30) 



which is in perfect agreement with the standard non-relativistic framework (see 
e.g. Ref. [29]). Here \x denotes the scale of dimensional regularization. The 
divergence at D — > 4 is ultraviolet. However, as mentioned above, in the following 
we shall not distinguish between the ultraviolet and infrared poles. Note that, 
since this expansion does not move the location of the singularity in the self- 
energy fixed by the relativistic relation p 2 = M 2 , there is no need to eventually 
re-sum the obtained series. 

The expression for the Z-factor of the charged pion can be directly obtained 
from Eq. (28) by performing the limit f2 — > at d > 3 (see, e.g., Refs. [30,35]). 
In this manner, one finds 

Z=l. (31) 

4.3 Coulomb vertex 

The Coulomb vertex, shown in Fig. 3b, is given by 

T/ 2 f d ° l Q ( ,j\q/ n (Wl + W[)(W 2 + W' 2 ) 

Vc = ~ e J 5 ±Cpi + *)S±(P2 - p , 

wi = w±(pi), w[ = w±(pi + 1) , w 2 = w±(p 2 ), U>2 = W ±(P2 - 1) • (32) 

Performing the Cauchy integration over 1° and using the threshold expansion, 
one can demonstrate that the numerator may be replaced by 

( Wl + w[){w 2 + w' 2 ) -> 4w lW2 - ( Wl - W2 Q + 2q\ 1 - l 2 , (33) 



3 Note that here and in the following, the results we show for various diagrams are always 
to be understood with the threshold expansion applied. 
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with Q = pi + P2 and q = (pi — P2)/2. Further, it can be checked that the 
energy denominator that emerges after integration over 1° can be rewritten as 



+ 



2w[2w 2 w[ + w' 2 -E 2E\ 2 + 2bl - (Ql) 2 /£ 2 - \(E 2 , w 2 , w 2 )/4E 

t? u Q(Qq) ,a A \ 

E = w 1 +w 2 , b = q ^— . (34) 

Using the substitution 

Q(Q1) Q(Qb) 2 2 2 

l-,l + -^ T -^-b-— ^ , s = £ - Q = (pi +p 2 ) , (35) 

we may rewrite the integral as 

v - e2 I dH 1 4M - + m 2 m\ 

Vc --%fij (2^P-,5(l-,„q/|q|) 2 + "" *-4- M " (36 > 

where the ellipses again represent higher-order terms in e. At lowest order in e 
we obtain (cf. e.g. Refs. [29,30]) 

Vfc=^(-^-Wc(«b)+0(d-3))+..., g 2 >0, (37) 
where a denotes the fine structure constant, 

^»^( ifa2A+ K in f (38) 

is the (infrared- divergent) Coulomb phase, and the ellipses stand for higher-order 
terms in the expansion in e. Since these terms do not change the position of the 
branch point in the diagram at s = 4M 2 , there is no need to re-sum these terms. 

From Eq. (36) we may observe that the counting of the photon momentum 
1 = 0(e) is different from the case of the self-energy diagram. 



4.4 Transverse vertex 

The vertex with the transverse photon, which is shown in Fig. 3c, is suppressed 
by a factor e 2 as compared to the Coulomb vertex due to the derivative coupling 
of the transverse photons to pions. It is given by the following Feynman integral 

Below, instead of Vp, we shall present the result for V c + Vr- Applying threshold 
expansion, it can be shown that the sum of Coulomb and transverse vertices can 
be rewritten as 

Vc + V T = e 2 J -£L S ± ( Pl + l)S ± ( P2 - I) ^l±l. . (40) 



16 



Using the same technique as in the case of the Coulomb vertex, we finally obtain 
a compact expression valid to all orders in the e-expansion, 



V C + V T = - - — + iO c {q Q ) + —= + G(d - 3 

= ^{^) +j (^ Waj) _0} +O(d _ 3 ), (4!) 

where in the second line, we have introduced a = i/l — AM^/s. As in the 
Coulomb vertex, the photon 3-momentum in the loop 1 counts as 0(e). Note 
that Vc + Vt is explicitly Lorentz-covariant. Equation (41) can also be derived 
in the formalism of Refs. [32,33], see Ref. [36]. 



4.5 Transverse vertex in crossed channel 

The vertex in which one of the photon lines is attached to the incoming kaon line 
is depicted in Fig. 3d. Neglecting all factors in the numerator, we consider the 
following scalar integral 

Vk = e2 / s±iPl - 1)Sk{pk ~ l) ^p- (42) 

After integrating over 1° we find 

K ~ 6 J J2n)~ d 2|T| ( w ' 1 - Wi + \\\)( w ' k - Wk + \\\) ' { A) 

where w[ = w±(pi — 1), w\ = w±(pi), w' K = wk(pk — 1) and wk = wk(pk)- 
From Eq. (43) one concludes that the photon momentum 1 has to be counted as 
0(e 2 ). However, since 

w' 1 -w 1 + \l\ = |1|-HL + ... u ,^_ u , A:+ |i| = |i|_P^d + .. (44) 

Wi W K 

and since the second term in this expansion counts as 0(e 3 ), we arrive at no- 
scale integrals. Consequently, we conclude that V K = 0; see also the discussion 
in Appendix E. 



4.6 Coulomb exchange in strong loops 

The exchange of a Coulomb photon in the strong loop (see Fig. 3e) is described 
by the following integral: 



Jc = e 7 s±{l)s±{Q ~ l)s ^ k)SM 

(45) 
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where Q = pi +P2, w[ = w±(l), w' 2 = w±(Q — 1), w'{ = w±(k), and w'{ = w±(Q — 
k). Replacing (w[ + w'{){w' 2 + w'^) = 2w±w 2 + 0(e 2 ), where W\ = w±(pi) and 
W2 = w±(p2), integrating over Z°, k° and re-defining the integration 3-momenta 
1, k along similar lines as in Sect. 4.3, we finally find 

_ e 2 M 2 C d d l d d k 1 1 1 

° ~ ~r J (2^p (2^p k^f + • " ' { > 

where q 2 = s/4 — M 2 and s = Q 2 . From the above expression it is seen that 
1, k count as 0(e). The higher-order corrections in Eq. (46) do not change the 
location of the threshold given by the relativistically invariant relation s = AM 2 . 
Calculating the integral in Eq. (46), we arrive at 

*~£(»^+K-2))+-- <«> 

where 

u 2(D ~ 4) / 1 \ 



5 Power counting 

We start by summarizing the power counting rules for K — > 37r decays without 
photons. 

1. Polynomial interaction terms can be organized in even powers of e, start- 
ing at e°. The inclusion of higher-order derivative terms in loop diagrams 
increases the power of e of the respective graphs accordingly. 

2. In loop diagrams, each pion propagator counts as e~ 2 , each loop integration 
as e 5 (as for an integration momentum /, d\ = 0(e) and dl° = 0(e 2 )). 

3. Each loop formed with two-body rescattering terms (counted as 0(a)) there- 
fore increases the power of e by e 5 (e~ 2 ) 2 = e 1 , therefore one-loop diagrams 
are 0(a 1 e 1 ), two-loop diagrams 0(a 2 e 2 ), etc. Loops with three-body terms 
increase the power of e by (e 5 ) 2 (e~ 2 ) 3 = e 4 . 

The examples of the preceding section show how these rules have to be amended 
in the presence of (virtual) photons. Note that, as M 2 — M 2 = 0(e 2 ) (when 
neglecting a tiny contribution of second order in the light quark mass difference), 
one might infer that the electric coupling should be included in a unified power 
counting, where e = 0(e). However, we do not follow this approach, but instead 
count e as a separate expansion parameter, such that we will altogether have 
a three-fold expansion in e, a, and e. In the following, we have to distinguish 
between "soft" photons (l = 0(e 2 ), 1 = 0(e)) and "ultrasoft" photons (l = 
0(e 2 ),l = 0(e 2 )). 
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4. Coulomb photons are always soft, hence their propagators are 0(e~ 2 ), and 
loop integration still counts as e 5 . Furthermore, the couplings to bosons are 
of O(e ). The addition of a Coulomb photon to a "skeleton" diagram, con- 
sisting only of interaction vertices without electromagnetism, hence changes 
the non-relativistic power of the diagram by a factor of e 5 (e -2 ) 3 = e _1 . The 
two most relevant examples for K — > 3n decays discussed individually above 
are the vertex correction with a Coulomb photon, Sect. 4.3, and Coulomb 
exchange within charged pion loops, Sect. 4.6. 

5. Soft transverse photons have couplings of 0(e l ), hence the exchange of trans- 
verse photons is suppressed by two orders with respect to the corresponding 
Coulomb exchange diagram. An example for this is given in Sect. 4.4. 

6. For ultrasoft transverse photons, with all momentum components of order e 2 , 
the propagator counts as e~ 4 and loop integration as e 8 . Therefore, adding an 
ultrasoft transverse photon to a skeleton diagram changes its power counting 
by a factor of (taking care of loop integration, boson propagators, photon 
propagator, and photon vertices) e 8 (e~ 2 ) 2 e~ 4 e 2 = e 2 . Examples are the pion 
self-energy (Sect. 4.2) and the exchange of a transverse photon in the crossed 
channel (Sect. 4.5). 

As the addition of a Coulomb photon modifies the power counting of a given 
diagram by a factor of e 2 /e, it is obvious that we can consider very small momenta 
for which this ratio ceases to be a small parameter. This is precisely what happens 
in the energy region where pionium bound states are formed and will be discussed 
in detail in Sect. 6. The scheme formulated above holds outside a small energy 
window around the pionium singularities, where then the standard perturbative 
treatment of electromagnetism is valid. 

As the inclusion of real photon radiation effects takes place on the level of 
widths or decay spectra, we briefly discuss the non-relativistic power counting 
for both non-radiative and radiative widths/decay spectra. The non-radiative 
differential decay width dT for K — > 3n is schematically given by 

dr oc [nwl^^-E^NI 2 - MPt) = 2p§^f- (49) 

Each invariant phase space element counts according to dfj,(pi) oc e 3 , while the 
energy- and momentum-conserving ^-function induces a factor of e -5 . With the 
matrix element Ai starting at O(e ), it is therefore obvious that the non-radiative 
differential decay width starts out as C(e 4 ), and the decay spectrum with respect 
to S3 obeys 

£ = W ■ (50) 
The radiative differential decay width gT 7 is schematically given by 

dT, oc \^d^ v S\d^k)5^( K P-Y J V l -k)\M 1 \\ dfi(k) = 2J f* n)3 ■ (51) 
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Real photons necessarily are ultrasoft, as the energy component is A; = 0(M K — 
3M n ) = 0(e 2 ), and the on-shell condition enforces the same order for the three- 
momentum components. Therefore the additional integration over the photon 
phase space is d/i(k) = 0(e A ). The non-relativistic power counting for the radia- 
tive matrix element is most easily seen in the Coulomb gauge, where the matrix 
element for photon radiation (with momentum k) off the external legs (with 
charges Q iy masses M iy and momenta p^ is given by 

M -r « e E^^( 1 + ( e )) = (£) • ( 52 ) 

Combining Eqs. (51), (52), we obtain dT 1 = (9(e 2 e 6 ) and 

dT 



7 = 0(e 2 e 4 ) . (53) 



We find therefore that, for the projected accuracy, it is sufficient to calculate 
bremsstrahlung effects at leading non-vanishing order in the non-relativistic ex- 
pansion. 

The following conclusions can furthermore be drawn. 

1. As the relation between the (massive boson) propagators in the relativistic 
and the non-relativistic theories is given by 

1 ' + O(e ) , 



p 2 -M 2 2co(p)(p° -w(p)) 

a calculation of the real photon radiation effects in both theories is equivalent 
up to the desired order. We can therefore simply expand the relativistic 
results in e. 

2. Contributions to the radiative width that are not of bremsstrahlung type, i.e. 
the photon is not attached to one of the external legs, is further suppressed 
in the e-expansion and does not need to be considered. 

3. The same goes for bremsstrahlung type contributions with derivative vertices 
or loops, which are also of subleading order in the e-expansion. 

At the end of this section, we can summarize our conclusions from the above 
considerations by discussing the Feynman graphs for all real and virtual photon 
corrections necessary for a consistent calculation of dT/ds^ up-to-and-including 
O(e 2 a e 4 ) (all channels) and 0(e 2 a l e 2 ) (for the "main modes" only). 

Virtual photon corrections suppressed with respect to leading tree terms by 
factors up to (9(e 2 e 2 ), i.e. without any nn rescattering contributions, are shown 
in Fig. 4. According to the counting rules spelt out above, the transverse photon 
diagrams Fig. 4a-c are all of the maximal order we consider here, so no further 
derivative vertices have to be taken into account for the same diagram topologies. 
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a) 




b) 




d) 




f) 




9) 



Figure 4: Representative set of one-loop diagrams for K — > 3ir including a virtual 
photon. Double, single, wiggly, and curly lines denote kaons, pions, Coulomb, 
and transverse photons, respectively. The thick blob denotes an interaction vertex 
with two derivatives. 



Figure 5: Diagrams for the bremsstrahlung process K — > 3^7. The line style is 
as in Fig. 4. 

We need to consider both diagram Fig. 4d for Coulomb (0(e 2 e~ 1 )) and Fig. 4e 
for transverse (0(e 2 e 1 )) photons. Coulomb exchange also necessitates the calcu- 
lation of the same graph with a two-derivative vertex, Fig. 4f, which is then also 
of 0(e 2 e 1 ). We furthermore wish to point out that no diagrams with photons 
coupling to interaction vertices, see Fig. 4g for an example, need to be consid- 
ered, as those necessarily contain transverse photons with derivative interactions 
and therefore only start to contribute at C(e 2 e 4 ). 

The bremsstrahlung diagrams are displayed in Fig. 5. As argued above, only 
their leading order in e contributes to dT/ds^ at (9(e 2 e 4 ), therefore no diagrams 
with additional derivative vertices are needed. In particular, structure-dependent 
photon radiation as in the graph Fig. 5c is suppressed by the absence of a meson 
propagator and the derivative in the vertex by at least two orders in e. 

Finally, we turn to the diagrams that contain one further tot interaction 
vertex, i.e. that are of O(e 2 a l e ). These are displayed in Fig. 6. Note that 
we only consider these corrections for the two "main modes" Ki — > 37r° and 
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a) b) c) d) 



Figure 6: Two- loop diagrams for K — > 37r, including a virtual photon and one 
mi rescattering vertex. Line style as in Fig. 4. 

K + — > 7r°7r°7r + . Obviously, the diagram Fig. 6a (discussed extensively in Sect. 4.6 
and below in Sect. 6) is the only one contributing to the former, and neither of 
the "main modes" features the topology shown in Fig. 6d (which requires at 
least two charged pions in the final state). As the pion loop as such is of 0(e), 
only the addition of a Coulomb photon can lead to the required order in the 
non-relativistic counting, and transverse photon contributions like Fig. 6c are of 
higher order and can be neglected. Finally, the graph shown in Fig. 6b could po- 
tentially contribute to K + —>■ tt 7t tt + , with a 7r + 7r _ pair in the loop rescattering 
into 7r°7r°. However, it is easily seen that the two contributions with the photon 
hooked to the tt + and the tt~ in the loop exactly cancel, hence we do not have to 
consider this diagram either. We conclude that the graph discussed in Sect. 4.6 
is indeed the only one at C(e 2 a 1 e°) that needs to be taken into account. 

6 Pionium 

In previous sections, we have analyzed the effect of radiative corrections in the 
perturbative regime, where an expansion in the fine structure constant a makes 
sense. In that region, pion momenta are counted as order e. Here, we investigate 
a region in momentum space where this is not necessarily true anymore: very near 
the cusp, pionium can form, and an expansion in a is obviously not possible. The 
three-momenta of pions are of order a, and one needs to take into account an 
infinite number of graphs to describe this momentum regime properly. 

There are two possibilities to cope with this problem. Either, one attempts 
to set up a formalism that allows one to treat the whole decay region, including 
the formation of pionium, in a coherent and consistent manner. Or, one proceeds 
in the manner in which the data analysis was performed by the NA48/2 collabo- 
ration [9]: one excludes a region around the cusp where this phenomenon occurs, 
and restricts the analysis to the region where a perturbative expansion in a is 
possible. We adopt this route in this article, for the following reason. Even if one 
might be able to set up a coherent formalism for the whole decay region, 4 there 

4 Note that, at present, there is no consent with respect to the pionium contribution in 
K — > 37r decays. For example, it was shown that, including this contribution, a better fit 
to the data can be achieved in the vicinity of the cusp [9]. However, the best fit value of 
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would be no real benefit of such a tour de force in the present context: as we are 
interested in a precise extraction of the irn scattering lengths, it is sufficient to 
stay in a region in momentum space where a perturbative treatment is possible, 
and where the effect of the cusp still allows one to extract the scattering lengths 
with high precision. Investigating the region where pionium is formed would 
then not provide any new information on the scattering lengths, but merely test 
the formalism to describe that region in a proper manner. This is so, because 
the energy resolution in present investigations of K — > 3tt decays does by far 
not suffice to e.g. measure the width of the ground state of pionium, which is 
an alternative possibility to measure nn scattering lengths. In our opinion, this 
alternative way is best covered by dedicated measurements of the production and 
decay of pionium as performed e.g. in the DIRAC experiment [18, 19]. Here, the 
theoretical background for the data analysis has already been provided with the 
necessary precision [40-43]. 

In the data analysis of NA48/2 [9], a region around the cusp was omitted, 

s t [AMI ~ A a , AMI + A a ] ; A s = 5.25 x 1CT 4 GeV 2 , (54) 

whenever corresponds to a channel where pionium can be formed. It therefore 
remains to answer the following question: is the region outside this interval ac- 
cessible to a perturbative expansion in a? As Ref. [9] has already demonstrated 
that this region does allow a precise measurement of the 7T7t scattering length 
(in the absence of radiative corrections), one is done once one knows that this 
question can be answered in a positive sense. The present section is devoted to 
an investigation of this problem. 

6.1 K L 3tt° 

It is useful to start with the decay K L — > 3ir°, because, as was shown at the end 
of the previous section, the graph Fig. 6a is the only one that needs to be taken 
into account in the perturbative region at order O(e 2 a 1 e ). It was evaluated in 
the present framework in Sect. 4.6. In the non-perturbative region, at order a, 
there are additional Coulomb exchanges that must be taken into account, see 
Fig. 7 for an example. The effect of these Coulomb exchanges is to modify the 
one-loop integrals J+-(s) displayed in Eq. (10), in the following manner, 

J+_(s) -> J + _(s) + J c (s) + Jnc(s) ■ (55) 

the K ± — > 71"^ + pionium decay rate normalized to the K ± — > n^Tr+ir^ decay rate is equal to 
(1.61±0.66) x 10 -5 and is considerably higher than the theoretical prediction <~ 0.8 x 10 -5 [9, 37] . 
Furthermore, it has recently been argued [38, 39] that this value substantially increases if the 
interference effect with the 7r°7r° intermediate states is taken into account, but the correction 
is model-dependent. A complete analysis of the problem is still pending [38, 39]. 
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Figure 7: Diagram with four Coulomb photons exchanged between the oppositely 
charged pions. The quantity J n c in Eq. (57) describes the sum of diagrams with 
n > 2 Coulomb photon exchanges. 

Here, the second term on the right-hand side denotes the contribution of Fig. 6a. 
At the order considered here, it is given by 

Ms) = (l6**A + In ]q l = \- Ml . (56) 

The third term on the right-hand side of Eq. (55) stands for the sum of n > 2 
Coulomb photon exchanges [44], 

JnC 

v 

Two remarks are in order. First, the ultraviolet contribution that shows up in 
Jc(s) can be absorbed in the couplings at tree level. In the following, we set the 
scale at \i — M n , and use instead of Jc the expression 

327T V M nJ 

The choice of the scale is irrelevant: a change of scale simply induces a correspond- 
ing change in the tree-level couplings. Second, we performed the evaluation of 
J n c in the non-relativistic framework described in Ref. [29, 30], which differs from 
the one proposed here by use of propagators where the non-relativistic expansion 
w±(p) — > M-x + p 2 /2M 7r + ■ • ■ is performed. A calculation with the propagators 
used here would be rather demanding. On the other hand, that calculation would 
lead to the same result at leading order in the momentum expansion. As we shall 
drop the contributions from multi-Coulomb exchanges anyhow, see below, we do 
not attempt to perform this considerably more complicated calculation. 

The effects of the non-perturbative regime is clearly seen in the multi-Coulomb 
contribution J n c- Indeed, once the pion momenta count as order a, v is of order 



(*(1 



*(!)) 



167T 



(57) 
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1, and a perturbative expansion is no longer possible. A result of this are the 
pionium poles that show up in J nC , at 



% 2 = -^-r ; n=l,2,.... (59) 

On the other hand, this phenomenon only happens for pion momenta that are 
sufficiently small. Expanding ^(1 — v) around v — 0, it is seen that multi- 
Coulomb exchanges are certainly negligible in comparison to the leading term Jc 
in case that the width of the cut satisfies A s > 10~ 4 GeV 2 . As a result of this, 
Coulomb exchanges finally amount to the replacement 

J+-(s) ^ J+-c(s) = J+-(s) + J c (s). (60) 

In the following, only the combination J + _c(s) will appear. 

We now investigate the size of the one-Coulomb exchange at order a and a 2 . 
We drop contributions from derivative interactions in the matrix element and 
evaluate the ratio 

dr int _ dr° 

R= dS3 ^ • (61) 

d,S3 

The quantity J^ mt is calculated with the replacement Eq. (60). We first consider 
the case where this replacement is performed at one-loop order only. Charged pion 
loops then generate the two graphs displayed in Fig. 8. On the other hand, the 
Coulomb contributions Jc are dropped in . The decay spectrum is calculated 
by use of the procedure described in Appendix F. 5 The result for the ratio R is 
displayed in Fig. 9, with a solid line, with tree-level couplings L = — \AK Q , see 
Eq. (8). The grey band denotes the region discarded in the NA48/2 analysis. It 
is seen that outside this region, the correction due to one-Coulomb exchanges is 
of the order of one percent or less. Note that the ratio R is linear in the ratio 
Lq/Kq to very high accuracy. 

We may also investigate R closer to the position of the pionium poles. In 
Fig. 10 we display the ratio R again, now also inside the region cut in the NA48/2 
analysis. The grey vertical band corresponds to S3 (jL [4M 2 — A s /5, 4M 2 + A s /5], 
i.e., the excluded region Eq. (54) is reduced by a factor 5. The dashed vertical lines 
correspond to the original cut Eq. (54). It is seen that also in this enlarged region 
of phase space, the electromagnetic corrections due to one-Coulomb exchanges 
remain small. 

We now comment on effects at order 0(e 2 a 2 e) - a contribution at this order 
is e.g. generated by the graph displayed in Fig. 11. A full calculation of all 



5 We note that, by using this procedure, the number of decay events in the vicinity of the 
cusp is increased by a factor three with respect to the standard evaluation of dT/d,S3. 



25 



+ 



Figure 8: A contribution to the ratio R in Eq. (61), evaluated with the replace- 
ment Eq. (60) at one-loop order. The pertinent R is displayed in Fig. 9 with a 
solid line. 




s 3 [GeV 2 ] 

Figure 9: The effect of one-Coulomb exchanges in the spectrum of Kl — > 37r° 
decays. Shown is the ratio R from Eq. (61). The grey band denotes the region 
Eq. (54), excluded in the data analysis of Ref. [9]. The solid line is evaluated 
with the replacement Eq. (60) at one-loop order. The dashed line corresponds to 
the same ratio, with the replacement Eq. (60) in bubble graphs at order a and a 2 . 
The dashed-dotted line illustrates that the effect of Coulomb exchanges at order 
aa 2 can be removed to a large extent through a renormalization of the coupling 
K in Eq. (8). 
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Figure 10: The effect of one-Coulomb exchanges inside the region Eq. (54) cut out 
in the NA48 analysis. The solid line again shows the ratio R at one-loop order, 
discussed in the text. The dashed vertical lines indicate the region Eq. (54), 
whereas the hatched grey band denotes the same region, decreased by a factor 
of 5. 




Figure 11: Contributions to the ratio R at two- loop order. In graphs with charged 
pion loops, one additional Coulomb exchange is retained. The corresponding ratio 
R in Eq. (61) is displayed in Fig. 9, with a dashed line. 
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pertinent graphs would be very complicated. For example, one would have to 
consider Coulomb photon exchanges in the two- loop graph displayed in Fig. la 
- a tremendous calculation. In order to have at least an order-of-magnitude 
estimate of the effect of these contributions, we perform the replacement Eq. (60) 
in bubble graphs at one- and two-loop order. The result for the ratio R is shown 
in Fig. 9, with a dashed line. It is seen that the change in R is rather substantial. 
However, almost all of this change can be absorbed in a redefinition of the coupling 
K in Eq. (8). This is illustrated with the dash-dotted line in the same figure, 
which results from the dashed one by changing Kq by half a percent, and which 
illustrates that the apparent sensitivity on the terms of order aa 2 is of no concern 
in this region of momentum space. 

6.2 K+ n°n°n + 

The discussion in this channel is analogous to the purely neutral channel just 
considered. We have checked that the corresponding ratio R - see Eq. (61) - 
has qualitatively the same behavior as is shown in Figs. 9 for the purely neutral 
channel, although with smaller absolute value. Furthermore, a renormalization 
of the tree-level coupling again allows one to remove most of the contributions at 
order aa 2 . 

We come to the following conclusions. The Coulomb corrections in hadronic 
loops are perturbative at least in the region which is used in the present data 
analysis by NA48/2 collaboration (7 energy bins discarded around the thresh- 
old). One may safely neglect pionium formation here. Moreover, taking into 
account the O(e 2 ae ) correction suffices: the two-loop correction can be effec- 
tively removed by a finite renormalization of K — > 3n tree-level couplings. For 
completeness, we nevertheless display in the following section the amplitudes 
with Coulomb corrections partially included also at two loops. According to the 
above discussion, we expect that the inclusion of these additional terms does 
not introduce a significant modification of the extracted values of the scattering 
lengths. 

7 Results: amplitudes and correction factors 

In this section, we collect the main results of our investigation: the correction 
factors Q and the modified matrix elements Ai mt for all four K — > 3tt channels. 
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7.1 K, 3tt° 

We start with the simplest case, with no external charged particles. Consequently, 
in the decomposition 



dT_ 

ds 3 



dT int 

^ooo -7— , (62) 
"S3 



where -E^ax i s ^he maximum photon energy in this channel, the external correction 
factor is just equal to one, 

^000 — 1 • 

dT mt is given as 

dT int = N d$ 3 \M^ \ 2 (63) 

(see Appendix A for the Lorentz-invariant phase space element N d& 3 ), where the 
matrix element M^g can be obtained from the matrix element Mooo in Ref. [7] 
by the replacement 

J + 4si) - J+-c(8i) , (64) 
see Eq. (60), which leads to the explicit form at two-loop order 

M$ = Mq cc + Ml- loop + A^^ loop , 

MT C = K + K 1 [(p° - M n o) 2 + (p° - M n o) 2 + (p° - M n o) 2 } , 

^l-loop = J 00 ( Sl ) + S (2)( Sl ) J + _ c ( Sl ) + ( Sl <- S2 ) + ( Sl <- S3 ) , 

_M^ loop = M£( Sl ) + M$( Sl ) + ( Sl <- s 2 ) + ( Sl <- s 3 ) , 
A<o (*i) = C 2 K J 2 ( Sl ) +AC X C+- L (J+-c(si)) 2 

+ (2C 2 x K + 2C 00 C x L )j + _ c ( Sl )J 00 ( Sl ), (65) 

and the polynomials in Ail~ loop as well as the genuine two-loop contri- 

butions M.q (si) are as given in Ref. [7]. Note that the replacement to J + _c*(sj) 
inside the two-loop contributions .Mjf (si) is in principle beyond the order of ac- 
curacy we are considering in this work, as it introduces a correction of 0(e 2 a 2 ), 
but it still represents a valid partial higher-order contribution. 

7.2 K L -> 77+77-77° 

In the decay channel Kl — > 7r + 7r~7r°, all radiative corrections are contained in 
the correction factor in the decomposition 

dT 



ds? 



dT int 

= ^ + -o(s 3 ,E-J—- (66) 



(where the expression for the maximum photon energy -E™ ax in this channel is 
as given in Eq. (A. 5)), as Ai mt determining <fT int only contains electromagnetic 
contributions of order e 2 a 1 , which we ignore in this channel. Ai mt is therefore 

Af + nt _ = M + - (67) 
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as given directly in Ref. [7]. Q + -o(s 3 , -E max ) comprises effects of virtual photon 
exchange and real photon radiation, which we may write as 

fi+-o(s 3! = 1 + £ (u™\(s 3 ) + u^ (s 3 , £ ma j) + 0{a 2 ) . (68) 

The virtual photon corrections also contain ultraviolet divergences, which are 
absorbed in a redefinition of the K — > 3ir coupling constants as defined in Eqs. (7), 
(8). In principle, these are of a form 

U = U + e 2 L t + 0(e 4 ) , U = L[(/i) + /3,A , (69) 

(and similarly for the Gi, Hi), i.e. the couplings Li (Gi, Hi) contain ultraviolet 
divergences and a scale-dependent part proportional to e 2 . We do not indicate 
these divergences explicitly, but remove all poles at D — 4 automatically. 

The Feynman graphs for the virtual-photon corrections are displayed in 
Fig. 4b, d-f. In the non-relativistic framework, we find the following result, 
using Eq. (41): 



ReA^+- = {^o + MPs - M - )} 



l + e 2 



a 



16(7 



(70) 



with a = a/1 — 4M2/s 3 . It is sufficient to determine the real part Re.M+_o, as 
the imaginary part would only contribute in |.M+-o| 2 by interference with the 
imaginary parts of ixix rescattering graphs and hence be of order e 2 a 1 , which is 
beyond the desired accuracy for this channel. Equation (70) translates into 

7r 2 (l + (7 2 ) 



virt 



o(53) 



(71) 



2(7 

Bremsstrahlung corrections in Kl(P) — > 7r + (pi)n~ {^2)^ \pz)l{k) are due to 
the diagram Fig. 5. The matrix element squared reads 

Ml Ml s, - 2Ml 



M: 



„2 7-2 

-e L 



+ 



(72) 



(pik) 2 (p 2 k) 2 (pik)(p 2 k) 

and we have shown in Sect. 5 that we can perform the calculation of the differ- 
ential decay width in its relativistic form. It can be written as 



dT 



Nd<5> 4 m: 



(73) 



where the precise form of the four-particle phase space element N d$4 is given 
in Appendix A. The calculation of u>^ (s3, -E max ) * s discussed in detail in Ap- 
pendix C.2, leading to 



real / „ rpn \ _2 
0l S 3, Anax) = ^ 



ln ^nax 



1 + 



5 - 2 In 



1 + a/I^ 



8M Kl (M Kl - M w o)M w oE^ 
Ao 



+ 0(e 4 ,a), 
(74) 



Performing the soft-photon approximation in the phase space integration amounts 
to setting 5 — > in Eq. (74), see Appendix C.3. 
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7.3 K+ 7r°7r°7r + 

For this channel, the decomposition of the radiative corrections effects is given 
by 

dr_ 

ds 3 



dT int 

Q o + (s 3 ,E max ) — — . (75) 
as 3 



We start by discussing the correction factor f2 00+ (s 3 , E max ), where E max is given 
in Eq. (A. 5). We write 

fioo+(s 3 , E max ) = 1 + ^ (u™l(s 3 ) + wS£(s 3 , ^max)) + 0(a 2 ) , (76) 

in analogy to what was done for f2 + -o- 

The Feynman graphs for the virtual-photon corrections are displayed in 
Fig. 4a-c. In the non-relativistic framework however, as discussed in Sects. 4.2, 
4.5 and Appendix E, all these virtual photon contributions vanish in dimensional 
regularization after applying threshold expansion. Hence we find 

Os)=0, (77) 

and the correction factor f2 00+ (s 3 , E max ) is given exclusively in terms of real pho- 
ton radiation effects. 

The radiative decay K + (P) — > ti q (j>\)ti {p2)^ + {p^)l{k) is given in terms of 
the squared matrix element 



lA/f 7 \ - -e 2 G 2 



M' K Mi 2Pp 3 



+ 



{Pkf ( P3 ky ( P3 k)(Pk) 



(78) 



Following Appendix C.l for the calculation of the radiative decay spectrum, one 
obtains 



real 
^00+ 



' s Enmx ) = \ c /in - 3 + -(4 - 5 C ) Jl^J c 

,3, max; 6M |_ M 2| ^ C ^ V 

8M K (M K - M w )M w E max 
6 C = . (79) 

Performing the soft-photon approximation in the phase space integration amounts 
to setting 5 C -> in Eq. (79). 

The "internal" photon corrections in are again comprised in the re- 

placement 

J + _( Sl ) -> J+- C (si) (80) 
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Figure 12: The ratio R tot defined in Eq. (82), see text for more details. The grey 
band denotes the region Eq. (54), excluded in the data analysis of Ref. [9]. It is 
seen that the radiative corrections are very small in this case. 



(compare Eq. (60)) in the amplitude Aioo+ displayed in Refs. [6,31,45], which 
leads to the following explicit form at two-loop order, 



A/f int 



M% ee + M^ loop + M 2 N loop 



MT = G + Gx(p° - M.) + G 2 {pl - M w ) 2 + GM - plf 



M 



1— loop 
N 



M 2 N ~ loop 



B N1 (s 3 )J + - C (s 3 ) + .6^2(53)^00(53) + {-B^ 3 (si) J +0 (si) + (si <-> s 2 )} , 
Mi(s u s 2 , s 3 ) + Mn(s u s 2 , s 3 ) , 

M B N = 4H C x C + ^(J + „c(s3)) 2 + 2[G C 2 X + H C X C 00 ] J+-c(s 3 ) Ms 3 ) 

+ G C 2 J 2 (s 3 ) + {AG Cl J 2 (si) + (si <-> s 2 )} , (81) 

and the polynomials -Bjvi, -8^2, -Bjv3 in -Mjy -10015 as well as the genuine two-loop 
contributions Aif^ are as given in Refs. [6,31,45]. 

We conclude this subsection with the following comment. The correction 
factor Q o+ turns out to be close to unity, and dT mt /ds 3 ~ dT/ds 3 . To illustrate, 
we define 



R 



tot 



o dr int 
iZ oo+^7 



d,S3 



dS3 



(82) 



Again, let us perform the replacement Eq. (60) only in diagrams with no deriva- 
tive couplings, and at one-loop order. We display the resulting R tot in Fig. 12 in 
the region defined in Eq. (54). We conclude from this figure that the radiative 
corrections to the amplitude are small in this channel. 



32 



7.4 On the soft-photon approximation 

We wish to briefly discuss the quality of the soft-photon approximation of the 
bremsstrahlung calculation in the corrections factors fi+_o and f2oo+ as discussed 
in the previous two subsections. In Fig. 13, these are plotted as functions of S3 for 
a fixed photon cut energy E cut = 10 MeV. Full lines denote the complete calcu- 
lation, dashed lines the soft-photon approximation. In addition, the dotted lines 
show a "naive" soft-photon approximation with E max = E cut , i.e. one disregards 
the fact that the maximum photon energy is, in some parts of phase space, not 
given by E cut , but by the kinematical limit E^^ss), see Eq. (A. 5). In the other 
cases, the correction factors are not smooth at the point where E^ n (s^) = E cnt . 

For the channel K L — > 7r + 7r~7r°, we find that, as power counting suggests, 
the Coulomb pole term in cu™* (s 3 ) completely dominates fi+-o- ^ n the i nser t 
of the upper panel of Fig. 13, this term has therefore been subtracted to make 
the differences between the various curves more visible. Even so, we find that 
the soft-photon approximation is a very small effect, and the difference in f2+_o 
stays well below 2 x 10~ 3 . However, there is a more significant deviation from the 
"naive" soft-photon approximation due to a logarithmic divergence at the upper 
kinematical limit for S3, signalling that very close to the boundary of phase space, 
the approximation of 0(a) becomes unreliable. 

For the channel K + — > 7r°7r°7r+, the correction factor is very small throughout 
phase space (|1 — fioo+l < 4 x 1CT 3 ), and the deviation due to the soft-photon 
approximation is smaller than 0.2 x 10~ 3 . Furthermore, except for the region 
where the maximum photon energy is restricted by phase space, the correction 
factor is very close to linear and could therefore be largely absorbed in a coupling 
constant shift of the polynomial part of the amplitude. 

We therefore confirm for both these K — > 3n decay channels that a treat- 
ment of the bremsstrahlung corrections in the soft-photon approximation leads 
to negligible modifications in the pertinent correction factors. 

7.5 K+ TT+TT+TT- 

As the corrections to the soft-photon approximation for both K L — > 7r + 7r~7r° and 
K + — > 7T + 7r 7r° turn out to be negligible, we confine ourselves to this approxima- 
tion for the last channel K + — > 7r + 7r+7r~. In this case, the universal radiative 
corrections cannot be formulated as a function of the single variable S3 any more, 
therefore we here define f2 ++ _(si, s 2 , s 3 , E max ) in terms of the differential decay 
width instead of the decay spectrum, 

dT\ = fi++_(si, s 2 , s 3 , E max )dT int , (83) 

n++_(si, s 2 , s 3 , £ max ) = 1 + ^ s 2 , s 3 ) + u£g_(s 3 , £ max )) + 0{a 2 ) . 
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Figure 13: Comparison of the photon correction factors in the full calculation 
(solid lines) and the soft-photon approximation (dashed/dotted lines), for the 
channels Kl — > 7r + 7r~7r° (upper panel) and K + — > 7r°7r°7r + (lower panel), as 
functions of s 3 . We have employed a photon cutoff energy of E cut = 10 MeV. 
In the upper panel, the main plot shows f2 + _ (s3, -^max) ( m units of 1CT 3 ), while 
the insert represents fi^ (s 3 , i?™ ax ) = a/n x uj™^ (s 3 , E^) alone, i.e. with the 
Coulomb pole subtracted. The lower panel shows f2 o+(s3, -E max ) (in units of 
10~ 3 ). The vertical lines denote the allowed kinematical range in S3. For more 
details, see main text. 
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We remark briefly that the definition Eq. (83) (in contrast to Eqs. (66), (75)) 
is only unambiguous as long as the soft-photon approximation is applied, which 
neglects the photon in overall momentum conservation. Otherwise the allowed 
kinematic range for, say, s\ and s 3 is different for the radiative process as com- 
pared to the non-radiative one. 

All diagrams shown in Fig. 4a-f contribute for this channel, and we find the 
following result for virtual-photon corrections: 



ReA4 ++ _ = {tfo + #i(P3-^r)} 



1 + e 2 



16(71 



+ 



1 + af 

16(79 



l + aj 

16(73 



(84) 



where we use the notation oi = \J\ — AM%/si. This results in 



7T 

w™ + _(si, s 2 , s 3 ) = — 



+ 



&2 



1 + al 

^"3 



(85) 



The soft-photon radiation contributions are calculated from the squared matrix 
element 



++- 



_ 2 tt2 \^ n n PiPj 
e ^0 Z^^^^.fcVr,, 



i,j=0 



{Pik){pjk) ' 



(86) 



where we use po = P and the charges Qo = Q3 = —1,Q\ = Q2 = 1 for 
compactness. The significantly simpler calculation results in 



LO 



real 
++- 



(87) 



The term +0(E ma , x ) indicates once more that, in contrast to earlier sections, real 
photon emission is only treated in the soft-photon approximation. 

As for the other supplementary channel K L — > 7r + 7r - 7T°, we disregard radiative 
corrections of order e 2 a 1 , therefore the amplitude M 1 f + _ is given directly in terms 
of the amplitude M.++- in Ref. [6,31,45]. 



8 On the accuracy of a determination of ao — a<i 

In this section, we discuss the various sources of theoretical uncertainties that 
limit the accuracy of an experimental determination of ao — a 2 from an analysis 
of the cusp in K — > 3n. For definiteness, we concentrate on the charged kaon 
decay channel K + — > 7r°7r°7r + . The authors of Ref. [4] were the first to quote a 
theoretical uncertainty on ao — 02, attributed to higher-order (three- loop) terms in 
the expansion in a as well as radiative corrections. A simple order-of-magnitude 
estimate of terms of 0(a 3 ) led to a combined uncertainty of 5%. This assessment 
was roughly confirmed in Ref. [5], where the uncertainty is quoted as 5-7%. 
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The analytic structure of the various K — > 37r amplitudes is now rather well 
understood, and combining the results of Refs. [6, 7] with the present investiga- 
tion, one obtains a very accurate representation. We identify the following main 
sources of theoretical uncertainty, in partial agreement with Refs. [4, 5] . 

1. Electromagnetic corrections. 

2. Isospin-breaking corrections in the matching of the coupling constants Cj to 
the 7T7T scattering lengths. 

3. Terms of order a 2 e 4 . 

4. Higher-order loop corrections in the strong sector, starting at C(a 3 e 3 ) in the 
non-relativistic power counting. 

Concerning point 1, we believe that the results presented in this article leave no 
significant error due to missing radiative corrections. This clearly is the main 
progress compared to earlier investigations. Concerning point 2, we have shown 
in Sect. 3.2 that the uncertainties are of the order of one percent. Examples of 
terms at 0(a 2 e 4 ) mentioned in point 3 can be worked out by comparing the exact 
form of the two-loop integral F in Eq. (13) with the explicit expressions at order 
e 2 , provided in Ref. [7], or by partially adding effective range interaction vertices 
to the two-loop diagrams. 

Finally, for an estimate of the uncertainty in point 4, we consider the thresh- 
old theorem [2, 3, 6] for K + — > 7r°7r°7r + in the absence of photons: the leading 
non-analytic piece oc a at the 7r + 7r~ threshold is given by the product of the 
decay amplitude K + — > 7r + 7r + 7r~ and the scattering amplitude 7r + 7r~ — > 7r°7r°, 
both evaluated at threshold. In particular, the coefficient of the cusp oc a is 
proportional to 

Z = T(K + -> 7r + 7r + 7r-)| lf . t2 „, 2 , . (88) 

A calculation of the quantity Z to 0(a n ) therefore yields the strength of the 
leading cusp behavior as given in a representation of the decay amplitude K + — > 
7T 7r°7r + at 0(a n+1 ), in other words, the two- loop representation of the K + — > 
7r + 7r + 7r~ amplitude given in Ref. [6,31,45] allows for an estimate of a? effects on 
the cusp. Numerically, the expansion of Z up to two loops reads in arbitrary 
normalization 6 

Z = -1.0 (tree) - 0.13 i (1-loop) + 0.014 (2-loop) , (89) 

i.e. terms of order a 3 will modify the strength of the leading (one-loop) cusp by 
about 1.5%. The present estimate does not replace a full calculation of C(a 3 ) 
effects, as it is necessarily incomplete: it neither yields a correct representation of 

6 We use a rough estimate for the couplings Go, Ho worked out from the slopes provided 
in PDG [46] and drop all derivative terms proportional to the couplings G\ , H\ , . . . , which 
generate small corrections only. 
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the amplitude elsewhere in the Dalitz plot (except in the cusp region), nor does 
it give any information about subleading cusp behavior (e.g. oc a 3 ). However, we 
regard Eq. (89) as a good indication for the rate of convergence in the K + — > 
7r°7r°7r + amplitude and the error in the latter due to the omission of (strong) 
three-loop effects. 

We may conclude from the above that the decay amplitude for K + — > 7r°7r°7r + 
is indeed known very accurately. However, the accuracy to which the amplitude is 
known does not necessarily translate directly into the accuracy with which a — a 2 
can be extracted from data. Rather, the latter is related to the change in the 
derivative of the amplitude (or the decay spectrum) in the vicinity of the cusp. 
We illustrate this point with the help of the loop function J + _c(s 3 ), see Eq. (60), 
which is responsible for the "internal" radiative corrections in K + — > 7r°7r°7r + . 
Below 7r + 7r~ threshold, the singular part of it is given by 



1 AM 2 

(a + aloga), a=J - 1 . (90) 



167T V V S ' 3 

The derivative with respect to s 3 slightly below threshold is given by 



_d_ 

ds 3 



(a + a log a) 



16tt 



1 + » + (91) 



where higher-order terms in A/M 2 have been neglected. So while the term of 
0(a) may represent a very small correction on the absolute value of the ampli- 
tude, it changes the slope near the cusp, given by J + ~c 011 a smooth (polynomial) 
background, by a relative factor of 2aM 7T /\/ r K that diverges near threshold. For 
example, at the lower bound of the exclusion region in the NA48/2 experimental 
analysis [9], A = A s = 5.25 x 10~ 4 GeV 2 (see Sect. 6), this relative factor is as 
large as 9%. 

We do not doubt that a more reliable estimate of the theoretical uncertainty 
in ao — a2 is feasible by working it out along with the data analysis. As all the 
ingredients are available now, it would be useful to e.g. i) study the stability 
of a — a 2 with respect to changes in the choice of the exclusion region A s , as 
illustrated by the example above; ii) work out the influence of a 2 e 4 terms as 
described above; iii) investigate the effect of Coulomb corrections at one vs. two- 
loop order. 



9 Comparison with other work 

Isidori [25] has considered the universal soft-photon corrections in multi-body 
meson decays and applied this specifically to the decay channel K + — > n + n + n^ . 
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The results for the virtual photon corrections correspond to a relativistic calcu- 
lation, as it is discussed in Appendix D, i.e. they agree with our non-relativistic 
results up to polynomial terms (that can be absorbed in a redefinition of the 
polynomial part of the decay amplitude) and up to higher orders in the expan- 
sion in e. Note that these higher-order terms in e are not complete in the sense 
that there are non-universal corrections of the same order (some of which have 
been discussed as neglected higher-order diagrams in Sect. 5). The real photon 
radiation corrections are considered in the soft-photon approximation, the effects 
of which we have discussed in Sect. 7.4. Our results agree with those given in 
Ref. [25] (up to typos) when expanding once more in e. Finally, Isidori [25] men- 
tions the formation of Coulomb bound states (such as pionium) as one of the 
potentially relevant omissions in its discussion of radiative corrections. We have 
closed this gap here in Sect. 6. 

In recent papers [26,27] the isospin-breaking corrections to K ± decays into 
three pions have been studied in a combined approach, where the treatment of the 
cusps resembles that in the non-relativistic effective field theory (see Refs. [6, 7] 
and the present work), whereas the electromagnetic effects are included within 
a quantum-mechanical approach. The treatment is not systematic. For exam- 
ple, in the strong sector the expression for the decay amplitude in our language 
corresponds to the summation of the one-loop pion loops with non- derivative 
couplings only. This expression has not been derived in Refs. [26,27] and, as 
follows from the present work, can be valid only in the vicinity of the cusp and 
only up to and including 0(e). The generalization to higher orders and beyond 
the cusp region is not discussed. The electromagnetic effects include the multi- 
Coulomb exchange within the charged pion loop that corresponds roughly to our 
Eq. (55). The procedure is again not systematic: Bremsstrahlung corrections, in- 
frared divergences, etc. are not discussed. We therefore conclude that the results 
of Refs. [26,27] cannot be safely used for the analysis of the experimental data 
on K — > 3n decays. 

The calculations of radiative corrections to K — > 3ir decays in Refs. [22-24] are 
performed in the framework of Chiral Perturbation Theory. Such a framework 
is not suited to an extraction of the tttt scattering lengths from experimental 
data, as these do not appear directly as parameters of the theory. Furthermore, 
the effects corresponding to our corrections of 0(e 2 a) that modify the analytic 
structure of the cusp are beyond the accuracy considered in Refs. [22-24] and are 
not considered there. 

10 Summary and conclusions 

In this article, we have generalized the framework of non-relativistic effective 
field theory for the cusp analysis in K — > 3n decays to include real and virtual 
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photon effects. This framework now allows one to calculate the relevant ampli- 
tudes systematically in a three-fold expansion in the non-relativistic expansion 
parameter e, titt scattering lengths a, and the electromagnetic coupling e 2 . We 
have performed this calculation up to C(e 2 a°e 4 ) for all four K — > 3n channels, 
and up to 0(e 2 a 1 e 2 ) for those channels in which the cusp effect is seen, namely 
K + — > 7r°7r°7r + and K L — ► 3ir°. We have made consistency checks with the cor- 
responding relativistic calculations of radiative corrections in various ways. Real 
bremsstrahlung has been calculated without the soft-photon approximation for 
K + — > 7r°7r°7r + and Kl — > 7r + 7r~7r°, and the effect of this approximation could 
explicitly be shown to be small. The non-relativistic power counting provides a 
natural explanation for the common assumption that radiative corrections are 
dominated by Coulomb photon exchange, and that e.g. bremsstrahlung effects 
are very small. As an important phenomenon beyond the universal radiative cor- 
rections, we have studied the impact of pionium on the threshold (cusp) region. 

We have provided simple analytic formulae for universal multiplicative correc- 
tion factors on decay spectra or differential decay widths, and for modified decay 
amplitudes taking into account "internal" photon exchange processes. We have 
commented on the various issues concerning the theoretical uncertainties in the 
extraction of ao — a>2 from the cusp effect. We expect that using this formalism 
in a refined analysis of the experimental data will lead to a precise determination 
of the 7T7T scattering lengths from K — > 3ir decays. 
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A Notation, kinematics 



In this appendix, we briefly collect some of the notation that is used throughout 
the article. The masses of the charged pions are denoted by while for the 
neutral pion, M w o is used. Accordingly Mk denotes the charged kaon mass and 
Mk l that of the Kl- The decay channels are abbreviated by indicating the charge 
of the final state pions, and the momenta are assigned in the usual convention 
such that the odd pion carries four momentum p 3 . 



decay channel 


abbreviation 


Kl(P) 




(0,0,0) 


Kl{P) 


-> 7r + (pi)7r-(p 2 )7r°(P3) 


(+,-,o) 


K + {P) 


^n ( Pl )n°(p 2 )n+(p 3 ) 


(o,o,+) 


K+{P) 




(+,+,-) 



The Mandelstam variables are defined as 

Si = (P2+P3) 2 , cycl. 



(a.i; 



In the charged kaon rest frame, the energy and 3-momentum of the pions are 
given by 



o M\ + M 2 - s t 

Pi 



2M, 



K 



2 _ \(M K , M?, 
Pi AM 2 K 



(A.2) 



(in the neutral kaon decays Mk — > Mk l )- In the above equation, we have used 

X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz . (A. 3) 

Furthermore, the following functions of kinematical variables are frequently used: 
Ao = A (M 2 Kl ,s 3 , M 2 ) , A c = A (M|, s 3 , M 2 ) , 

(A.4) 



AM 2 



S3 



S3 cr 



The maximum photon energy S max (-^wJ f° r K + (Kl) decays in the kaon rest 
frame is given by 



E r 



(n) 
kin 



> Scut 



cut 



v cut • -^kin 

M 2 K - (M n + ^T 3 ) 2 



2M 



K 



rpn 

-^kin - 



where 

M 2 Kl - (M n o + ^~ 3 ) 2 



2M, 



k l 



(A.5) 
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Loop integrations are denoted by an angle bracket and the pole in four di- 
mensions is contained in the term A, 



(A.6) 



167T 2 



D-4 2 



i(r'(l) + ln(47r) + l) 



In the text, we also use - depending on the context - the symbols Xjr or X uv for 
this quantity. 

The differential decay width for the n-body decay of a particle of mass M in 
D = 4 dimensions is given in terms of the matrix element M. according to 

dT = Nd^ n (P; Pl ,...,p n )\M\ 2 , 



N 



2M 



■-, dMP- Pl ,..., Pn ) = 5^P-J2p^U^y- (A.7) 



where the values of the symmetry factor and the mass, (S,M), are (2, Mr-) for 
(0,0,+) and (+,+,-), (1, M Kl ) for (+,-,0) and (6, M K J for (0,0,0). For the 
phase space integrations inD = d+1^4: dimensions, we generalize Eq. (A.7) 
using 



AT/JS (2ir) D 1 J . . d d p 2tt^ 



2M S ' ^ 7 2p°(27r) rf ' a T(^±i) ' 

n n 

d$ n (P; Pl , . . . , p n , d) = 5 D (P - Pi) II • ( Ai 



i=i i=i 



The cut of the logarithm is chosen, as usual, on the negative real axis and the 
ip function and the dilogarithm are defined as 



Mx) = -^-]nF(x), U 2 (z) = - C dx ln{l XZ) . (A.9) 

dx J x 



B Radiative corrections: relativistic framework 

The non-relativistic calculation presented in the main text is quite complex. In 
order to have a fully independent check on that calculation, in this appendix 
we evaluate radiative corrections in a relativistic framework, for the two decay 
channels (+, — , 0) and (0, 0, +). In particular, we work out the corrections fi + _ 
and fioo+ ( see Eqs. (66), (75)) in the relativistic framework and compare the 
result with the non-relativistic calculation described in the main text. We confirm 
that, to the order considered in the momentum expansion, the relativistic result 
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indeed is identical to the non-relativistic one, up to a redefinition of tree-level 
couplings. Note that, as the purpose of this Appendix is mainly to cross-check 
the non-relativistic calculation, we also here disregard all one-photon and one- 
pion reducible graphs (see e.g. Ref. [23]), which turn into polynomial terms when 
expanded non-relativistically in e. 

In the present Appendix, the final results of this calculation are presented. 
The pertinent contributions include virtual and real photons. The evaluation 
of Bremsstrahlung is performed in Appendix C, whereas the relativistic loop 
integrals are detailed in Appendix D. 

The calculation is performed in the Feynman- gauge. Infrared and ultraviolet 
divergences are both tamed with dimensional regularization. In keeping with 
the standard procedure for calculations of radiative corrections in a relativistic 
framework, we retain the distinction between both types of divergences in this 
appendix. To make the formulae more transparent, we use the same scale for 
both, UV and IR divergences, although the origin of the divergences can be 
clearly identified. We therefore use the notation 

,D-4 



A - ^ 

AUV — 



16tt 2 



1 1 



D-4 2 



-(r'(l) + ln(47r) + l) 



= Xir. (B.l) 



The loop integrals A+, A , Tf, J ab may be found in Appendix D, and the Brems- 
strahlung integrals // are given in Appendix C, together with the symbols <5o )C - 
As already mentioned in the main text, we evaluate real photon emission beyond 
the soft-photon approximation. 

B.l K + 7r°7r°7r + 

We first consider the decay K + (P) — ► n° (pi)n (p 2 )n + (p 3 ) . The relativistic inter- 
action Lagrangian is 

C mt = ^(l + e^^WV + h.c.) , 

b = Q\ uv + b r (n) . (B.2) 

Here and in the following, we do not explicitly write the couplings to photons, 
because these are trivially constructed via minimal coupling. The virtual correc- 
tions displayed in Fig. B.l generate the contribution 

M virt (s 3 ) =G [Z- e 2 A + ( S3 )] ; s 3 = (P - p 3 ) 2 . 

The Z-factor is 

Z = 1 - 4e 2 (\ uv - X IR ) . (B.3) 

Since the amplitude depends only on s 3 , integration over phase space gives for 
the decay spectrum in d dimensions 

f- = N(d)(2n)% 2 (Mls 3 ,mid)<S> 2 (s 3 ,mimld)^ (B.4) 
as 3 
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Figure B.l: Virtual photon corrections to the K + — > 7r 7r 7r + decay. 



where $2 denotes the two-particle phase space in d dimensions, see Eq. (C.4), 
and where pf = mf. We add Bremsstrahlung contributions, worked out in Ap- 
pendix C. These cancel the infrared divergences generated by Ai vnt . Finally, we 
determine the relativistic analogon of the correction factor introduced in Eq. (75). 
Using the notation Qqq + — > £^ 0+ , we find 



ft£ + (s 3 , £ max ) = 1 + -[Fx + F 2 + F 3 ] + 0{a 2 ) , 

7T 

F 1 = A(L E T 1 K *(s 3 )-T 2 K *(s 3 )) , 



(B.5) 
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max ) ± 3 (■S3, -^n 
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F 3 = 8vr 2 6 r (/i) - 3 In — - 2L E + 



where 



L E = ]n[ 2E ™^ K \ ^ I 



m^-8. 2 J-( S3 ) 



M|- + M 2 - s 3 



'\ c I Ml 
To compare with fi 00+ , we expand Qq 0+ in e up to and including C(e 3 ) 



(B.6) 
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M, 



-2 In 



1 + 



(B.7) 



It is seen that differs from f2oo+ by a polynomial in the external momenta, 
which can be eliminated by a redefinition of the tree-level coupling constants. We 
conclude that the relativistic calculation confirms the non-relativistic result for 
f2oo+ presented in the main text. 
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Figure B.2: Virtual photon corrections to the — > 7r + 7r 7t° decay. 



B.2 Ki 



7T+7T 7T° 



Here we consider the decay Kl{P) — > 7i + (pi)n (P2)^°{P3)- The relativistic inter- 
action Lagrangian is 



Ant = L Q (1 + e 2 a)K L 7T°7ih , 
a = 6A[/y + a r (/i) . 

The graphs displayed in Fig. B.2 generate the amplitude 

M™\s 3 ) = L [Z + e 2 A (s 3 )] ;s 3 = (P- p 3 f 



(B.8) 



(B.9) 



The Z-factor was already given in Eq. (B.3). We use the analogue of Eq. (B.4), 
and combine the result with real photon emission considered in Appendix C. The 
relativistic correction factor fi+_ , the analogue of fi + _ i n Eq. (66), becomes 



a 



1 + - (/i + h + h) + 0(a 2 ) 
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where 
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(B.10) 



(B.ll) 
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Expanding the result in e up to and including 0(e 3 ), we obtain 
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2^/1 -5o- 2 ln( 
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+ 0(e 



+ 0(a 2 ) . 



(B.12) 



This expression agrees with f2+_o in Eqs. (68), (71), (74) up to a polynomial in 
the external momenta, which can be absorbed with a redefinition of the tree- 
level couplings. We conclude that the relativistic calculation confirms the non- 
relativistic result for fi + _ presented in the main text. 



C Bremsstrahlung beyond the soft-photon ap- 
proximation 

In this appendix, we evaluate the phase space integrals used in the calculations 
of the bremsstrahlung corrections to K — > 3n decays, with a kinematical cut E cut 
in the photon energy in the rest frame of the kaon. The evaluation of the decay 
spectra is, therefore, performed in this frame. 

Most of the notation is set in Appendix A. In order to tame infrared singular- 
ities, we perform the phase space integrations in d dimensions and indicate this 
with an additional argument d in the infinitesimal n-body phase space element 
d$ n . The differential decay width is defined as 

dT = N(d)d^(P; Pl ,p 2 ,q,k,d) \Mf , (C.l) 

where q and k are the momenta of the odd pion and the photon, respectively. In 
a first step, the four particle phase space is separated in a three particle phase 
space integration over a two particle phase space integration, 

J d$ 4 {P;p 1 ,p 2 ,q,k,d) = (2n) d J ds 3 d$ 2 {Q;Pi,P2,d)d<5> 3 {P;Q,q,k,d) , 
where Q 2 = s 3 . 



C.l K+ 

We illustrate the calculation with the simplest possible case (0, 0, +). The ampli- 
tude squared is given in Sect. 7.3. After performing the two particle phase space 
integration, the integration over Q and over the angle between the photon and 
the charged pion as well as their angular parts, one is left with the integrations 
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over the photon- and the 7r + -energy. The boundaries of these integrations are 
given by 

X = A c + 4kM K A + , Ai = s 3 + Mli -M 2 K + M K k , 

_ M 2 K - (M n + ^r z f _ J £ kin : £ kin < E cut 

-C'kin — ^77 , ^max ~ \ rp rp ^ it. ■ 



2M^ ' [ ^cut : -Skin > -^cut 

The integral reads 

-f- = -G 2 e 2 (27r) d N(d)^ 2 (s 3 ,M 2 ,M 2 ,d)[I 2 + I 3 -2h} 
as 3 



/ dkk d ~ 5 / dq fi(k,q°,d), (C.3) 



8(2tt) m 

where the two-particle phase space volume in d dimensions is given by 

'\(s,m\,m 2 2 ) 



$ 2 (s,m?,m^) = J d$ 2 (P;Pi,P2,d) 



2 ^\ ^ 



(47r)*r(f)V5' 

(C.4) 

with P 2 = s, p 2 = m 2 , and 



/ 2 (<Z°,fc,d) (,/vf - ^)'' \ f 3 (q°,k,d) = M 2 {q{) _ (r)2 

U(q°, k, d) = g o (^EZt! , g = ^2 _ M 2 

g u — qz V 

M£ + Ml - s 3 + 2(kq° - M K q° - A;M X ) 



2kq 

We carry out the calculation in quite some detail for the integral I 2 , the other 
integrals can be performed along the same lines. 
Introduce the variable u as q° = q° + ku 

gd-ggg-i f E ™ „ Id _ 4 r +(fc) 

2 8(2vr) M 
with the functions 



/ dfcA;^ 4 / dug(u,k) d - 3 , (C.5) 



_ M| + M 2 - g 3 - q° ± -q° 

q= 2M K ' q= 2M^ U±{k) = ^— 



z= q°-M K (l + u) + ku^ ^—^ 
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To isolate the infrared singular part we split up the function g in a part with 
photon momentum zero and a remainder, 



g(u,k) d - 3 = g(u,0) d - 3 + Ag(u,k), 
Ag(u,k) = g(u,k) d - 3 -g(u,0) d - 3 

'g(u,k) 



= (d-3)\n 



g(u,0) 



+ 0((d-3) 2 ). 



(C.7) 



Since the photon momentum dependent function Ag is proportional to k a with 
a > 0, this integral vanishes in the limit d — > 3 and we can drop it. The phase 
space integral has an infrared divergence proportional to (d — 3) _1 , therefore, in 
a next step, we expand the remaining function g(u, 0) d ~ 3 in the vicinity of d = 3 
and keep only terms up to and including 0{d — 3). Therefore, the original integral 
reads 
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8(2tt) m * 
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(C.8) 



Splitting up the functions G 2 {k) and H 2 {k) as well in a photon momentum de- 
pendent and independent part, the original integral can be written as 



h = + ^7 J 2( S 3, ^nax) + 0(d - 3) . 



(C.9) 



where 
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The part I 2 R contains the infrared singularity of I 2 as well as finite terms. There- 
fore, the part can be evaluated in d = 3. Applying the same procedure to the 
remaining integrals, they divide again into two parts, 
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(C.11) 
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and one finds 
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The remaining integrations over u and k can be done numerically in three dimen- 
sions, as they are free of infrared divergences. 

The bremsstrahlung decay spectrum can be simplified considerably by ex- 
panding the infrared finite part of the combination I 2 + I3 — 2/4 that appears in 
Eq. (C.3) in the non-relativistic parameter e. We find 



# 2 (0) + # 3 (0)-2# 4 (0) = 



x 3 J 2 



3QM A K Ml 



5 - 6 In 



M K M W 



0(e 5 ) , 



/■En 

Jo 



dk 



AG 2 (k) + AG 3 (k) - 2AG 4 (k) 



k 



^3/2 
9M 



5r = 



6/A / 

|^^4-(4-5 c )v / T 3 ^ + 31n 
8M K (M K - MJMkEwk 



1 + 



+ 0{e>) , 



A, 



(C.13) 



C.2 K, 



7T + 7T 7T° 



In the decay channel (+, — , 0), the four particle phase space is again separated in a 
three particle phase space integration over a two particle phase space integration. 
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However, because the amplitude squared indicated in Sect. 7.2 depends on the 
momenta of the two particle phase space integration, this two particle phase space 
does not factorize anymore, but produces the functions R\ and R 2 , 
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where R\ and R 2 are evaluated in a vicinity of d = 3 and the two particle phase 
space volume $2 is given in Eq. (C.4). The remaining three particle phase space 
integrals can be solved along the lines of the previous section, and one obtains 
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where 
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Since the kinematics in the channel 
integration boundary has to be changed to 

(M n0 + 
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Again, expanding the infrared finite part of I\ in e leads to substantial simplifi- 
cations, 
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Furthermore, we find for the combination appearing in Eq. (C.14) 



1 - a 2 



Ri 
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1+ (ln2- |)(d-3) 



+ C(e 4 ). (C.19) 



Note that A = 0(e 2 ) and a = 0(e). 



C.3 Soft-photon approximation 

Performing the soft-photon approximation in the four particle phase space inte- 
gration amounts, in the present case, to discard the photon momentum in the 
5-function, 

dT = N(d)d$ 3 (P; Pl ,p 2 , q, d)dfi(k) \M \ 2 . (C.20) 

Note that energy and momentum are not conserved anymore. Comparing the 
exact result calculated in the previous Sects. C.l, C.2 with the soft-photon ap- 
proximation shows that 

i) the soft-photon approximation reproduces I- R in a vicinity of d = 3 

ii) the difference between the two results is analytic in E max and of 0(E m£LX ). 
The soft-photon approximation if for the integral Ii can be written as 

/? = //*+5^//(as,O) + 0(d-3), i = l,...,4. (C.21) 

For a discussion of the numerical size of the error induced by the soft-photon 
approximation in the non-relativistic limit, see Sect. 7.4. 



D Relativistic loop integrals 

Here we provide the loop integrals needed in the relativistic calculation presented 
in Appendix B. In particular, we evaluate the vertex functions A and A + , 
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generated by the vertex graphs displayed in Fig. B.2 and Fig. B.l, respectively. 
In this appendix, pion and kaon momenta refer to charged particles throughout, 



p 2 = Ml P 2 = Ml 
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where s = (pi + P2) 2 , t = (P — p 3 ) 2 . The evaluation of the integrals J ab is 
standard. In Appendix B, we use the notation J ab (p 2 ) = J ab (p 2 ) — J ab (0). The 
vertex integral is 
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For equal masses M a = M b = M w , we find 
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E Infrared divergences 



Throughout this paper, we have used dimensional regularization to tame both 
ultraviolet and infrared divergences. In dimensional regularization, applying the 
threshold expansion, one sets so-called no-scale integrals to zero. While this 
simplifies calculations considerably, the identification of infrared singularities can 
be complicated. We shall demonstrate this explicitly on the basis of the example 
of the vertex diagram, given in Figs. 3b+c. 

Let us start from the relativistic case. The pertinent scalar integral is given 
by the function G(s) = G(s, M 2 , M%) discussed in Appendix D. As shown there, 
the integral is ultraviolet-finite. We can rewrite Eq. (D.2) according to 

ReG(s) = -— In ( ] X IR + ReG r {s) 

scr \ 1 — a J 

1 



Ml 



A /i? + 0(q 2 )+ReG"», (E.l) 



where q is the relative momentum of the pion pair in the final state, o = 
a/1 — 4M 2 I s an d the quantity ReG"*(s) is finite as D — > 4. The quantity Xjr 
is defined by Eq. (A. 6). In order to emphasize that the divergence at D — > 4 is 
infrared, the subscript "IR" is attached. 

Let us now consider the same diagram in the non-relativistic effective theory. 
To this end, one may replace the pion propagators 

111 . , 

+ . . . (E.2) 



Ml - p 2 2M W Mjr + p 2 /(2M 7r ) - p° 

(for simplicity, we use here the standard version of the non-relativistic EFT). 
Performing the contour integration over 1°, at lowest order in the momentum 
expansion we arrive at 

r NR ( , J_ [J^^f^ 1 1 A 

{ ) AM W J (2tt)^1 2 VP-2ql 2M^ |1| - ql/M w + 1 2 /(2M W ) / 1 ' ) 
The real part of the first term is finite as D — > 4. Consequently, 

G ™ (s ^ " in I f i.| + .v(2m,) + °< q2 » + flnite at D - 4 • < E - 4 » 

Performing the remaining integration, we exactly reproduce the infrared diver- 
gence in Eq. (E.l). However, if the integrand is first expanded in inverse powers 
of M n , each term in this expansion is a no-scale integral and vanishes in dimen- 
sional regularization. Thus, if threshold expansion is applied, the non-relativistic 
EFT fails to reproduce the infrared divergences of the relativistic theory. 
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In order to understand this apparent contradiction, note that only the first 
term in the threshold expansion, which is singular at the origin as |1| -3 , is infrared- 
divergent. Introducing a splitting of the integration interval, one finds 

i f r B r c \ d d i 1 

= M^l mi) JSL.G^3 <*" ~ Ad ' 3) + ^ " Bd ' 3) ) 

= -^(X IR -X uv ) + 0(d-3). (E.5) 

In the second line, we have taken d > 3 (d < 3) in the first (second) term. We 
conclude that, although the above expression is formally zero, one may identify 
parts of this expression with infrared and ultraviolet divergences. The infrared 
divergences, which are present in the relativistic theory, are reproduced in non- 
relativistic EFT. 

To summarize, using no-scale arguments in dimensional regularization changes 
the structure of infrared divergences in the non-relativistic EFT, since it amounts 
to setting expressions of the type Xjr — Xuv to zero. It is further seen that the 
change in the amplitude is a low-energy polynomial, because of the presence Xuv 
in the "dropped" term. Then, owing to the fact that the infrared divergences can- 
cel at the end, one may justify using the shortcut, based on threshold expansion 
and removing all divergences at D — > 4 by the counterterms in the Lagrangian. 



F Phase space in Kl —> 37T 

In case of the decay K L — > 3ir°, it is convenient to take advantage of the symmetry 
of the final state which contains 3 identical particles [47, 48]. One introduces polar 
coordinates for the kinetic energies Tj = — M w o according to 



T x = ^(l + rcos0) , 




see Fig. F.la for the definition of the variables r, 0. In Fig. F.lb, we show the 
boundary of the physical region 

(1 + x)r 2 + xr 3 cos30 = 1, (F.2) 
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Figure F.l: The Dalitz-plot for K L — > 3ir° decay. In Fig. a), we display the polar 
coordinates, with p = Qr/3. The physical region for the decay is bounded by 
the curve displayed in Eq. (F.2), shown in b) with a solid line. The three dotted 
lines denote the location of the cusps at Sj = 4M 2 . By a proper choice of the 
outgoing momenta of the pions, all events may be mapped into the region abc. 



where Q = M Kl — 3M w o is the energy release, and x = 2e/(2 — e) 2 , e = Q/M Kl . 
The differential decay width is 



vi 



dT = ^M^ L d ^ {P] Ph P2 ' P3) lM]2 • (R3) 

Because the matrix element Ai is a Lorentz scalar, one can immediately perform 
seven of the nine integrations over the pion momenta, without a need to know 
the matrix element M.. Depending on the variables chosen, one has 

y"(27r) 9 d$ 3 = J C 23 ds 2 ds 3 = J C r4> r dr d<\> = J C s4> ds 3 d(j> , (F.4) 

with Jacobi determinants 

r _ - 2 r - 2 ^ 2 r _ « 2 V3 \Ml L + 3Ml -3s 3 \ 

C23 -4Mf L ' Cr *-— Q ' C * ~ 2AMl L cos 2 (0 - f) • (K5) 

We note that the momenta of the three outgoing pions can always be chosen such 
that the decay products are located in the triangle abc in Fig. F.lb. This is a 
perfectly legitimate procedure, because the matrix element is totally symmetric in 
the variables si, s 2 and s 3 in this channel. In case that the integration is restricted 
to the triangle abc, the pertinent Jacobi-determinant must be multiplied with 3!, 
because there are six identical contributions to the decay probability. In Fig. F.2, 
we show the differential phase space in the standard case where the full phase 
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0.14 



s 3 [GeV 2 ] 

Figure F.2: The differential phase space d^s/ds^. Its shape depends on the region 
chosen. The dash-dotted line is the standard full phase space region. The solid 
line corresponds to the case where all events are mapped into the triangle abc, 
see Fig. F.lb). 



space is used, and in the case where the events are mapped into the triangle abc. 
Although the integrated phase space is of course identical in the two cases, their 
shape is very different. In the main text, we use the variables s 3 , and restrict 
the phase space to the triangle abc. Note that by this procedure, the number 
of neutral pion pairs located in the vicinity of the cusp region is increased by a 
factor of 3. 
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